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1.  INTRODUCTION 


This  report  is  the  final  report  for  a  research  study  conducted  at 
Physical  Sciences  Inc.  for  the  Air  Force  Rocket  Propulsion  Laboratory  at  Ed¬ 
wards  Air  Force  Base  California  under  Contract  F0461 1-80-C-0042  to  determine 
the  absorption  coefficient  of  water  vapor/propellant  mixtures  for  laser  pro¬ 
pulsion  and  laser  heated  thrusters. 

A  laser-heated  rocket-thruster  converts  laser  energy,  beamed  from  a 
remotely  stationed  laser,  into  thrust  by  using  the  laser  radiation  to  heat 
the  propellant.  This  technique  provides  a  distinct  advantage  over  conven¬ 
tional  chemical  rocketry  by  yielding  a  high  specific  impulse  (greater  than 
1000  s  vs.  approximately  500  s  for  chemical  rockets)  at  high  thrust  levels 
(potentially  greater  than  1000  lb) ;  however,  significant  technological  ad¬ 
vances  in  the  design  and  development  of  high  power  laser  facilities,  pointing 
and  tracking  systems,  and  collection  optics  will  be  required  before  laser- 
heated  thrusters  leave  the  laboratory.  Either  pulsed  or  continuous  wave  (CW) 
lasers  can  be  used  as  the  radiation  source.  Examples  of  each  type  of  rocket 
are  shown  in  Fig.  1.  To  obtain  the  highest  specific  impulse  for  a  given 
temperature,  hydrogen  is  the  propellant  of  choice. 1-4  A  10.6  micron  Cf>2 
laser  is  a  potential  candidate  as  the  source  of  the  laser  radiation. 

The  only  absorption  mechanism  for  pure  hydrogen  at  this  wavelength  is 
inverse  hremsstrahlung,  which  requires  the  presence  of  electrons.  Ionization 
of  the  hydrogen  and  the  production  of  electrons  is  Initiated  only  at  tempera¬ 
tures  between  7000  K  and  11,000  K,  depending  on  the  pressure.  Thus,  pure  hy¬ 
drogen  must  be  heated  to  absorb  the  laser  energy. 

In  a  pulsed  laser  propulsion  system,  this  heating  is  accomplished  by 
focusing  the  laser  beam  to  create  a  laser-induced  breakdown.  Hie  hot,  rapidly 
expanding  plasma  drives  a  laser-supported  detonation  (LSD)  wave  into  the 
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surrounding  cool  hydrogen.5'6  This  gas  is  heated  to  extremely  high  tempera¬ 
tures  around  20,000  K)  and  ionizes}  and  the  electrons  generated  behind  the  LSD 
wave  absorb  the  laser  energy  to  sustain  the  propagation  of  the  wave.  The 
propagation  of  the  wave  continues  until  the  intensity  of  the  laser  beam  drops 
below  the  level  required  to  sustain  the  wave  or  until  the  pulse  ends.  The 
heated  gas  then  expands  in  the  rocket  nozzle  creating  thrust.  Laboratory 
scale  proof  of  concept  experiments  at  10.6  u,  yielding  specific  impulses  over 
1000  s,  have  been  conducted  at  PSI  and  have  demonstrated  the  feasibility  of 
this  system.5'6 

CW  laser-heated  thrusters  with  pure  hydrogen  propellant  operate  in  a 
similar  manner.  In  this  thruster,  when  laser  radiation  is  focused  into  a 
plasma,  a  laser-supported  combustion  (LSC)  wave  is  created.  If  the  gas  flow 
velocity  into  the  thrust  chamber  is  matched  to  the  wave  velocity,  then  a  sta¬ 
tionary  "flame"  front  is  created.  Hydrogen  passing  through  the  flame  front  is 
heated  and  ionized.  'Hie  electrons  generated  by  the  ionization  absorb  by  in¬ 
verse  Bremsstrahlung,  and  sustain  the  wave.  One  problem  with  this  type  of 
thruster  is  that  there  may  be  severe  radiative  heating  of  the  chamber  walls  by 
the  stationary  20,000  K  plasma  neceeaary  to  absorb  the  laser  radiation.  Seed¬ 
ing  the  hydrogen  with  an  easily  ionizable  substance,  such  as  cesium,  has  been 
proposed  as  one  method  to  lower  the  operating  temperature  of  the  thruster  and 
therefore  reduce  the  heat  load  on  the  walls.  If  the  hydrogen  is  seeded  with  a 
molecule  that  absorbs  at  low  temperatures  then  the  system  is  simplified. 

First,  no  LSC  wave  need  be  sustained  to  heat  the  hydrogen,  allowing  greater 
flexibility  in  the  setting  of  gas  flow  rates.  Second,  and  perhaps  more  impor¬ 
tant,  the  operating  temperature  is  substantially  lower  and  thus  the  heat  load 
on  the  chamber  walls  is  reduced.  The  problem  is  to  find  a  low  molecular 
weight,  low  terqperature  absorber  that  will  not  deconpose  to  any  significant 
degree  below  3500  K  to  4500  K  the  temperature  range  of  a  thruster  operating 
with  specific  ingsulses  between  1000  s  and  2000  s.) 
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Water  and/or  ammonia  have  been  proposed  as  possible  candidates  for  ab¬ 
sorbers  at  10.6  u.  Water  Is  particularly  attractive  since  significant  disso¬ 
ciation  does  not  occur  until  the  tesgjerature  exceeds  3600  K.  In  this  program 
measurements  were  made  on  the  10.6  u  laser  transitions  P(16),  P(18),  P(20) 
found  in  high  power  CW  C02  lasers  to  determine  the  absorption  coefficient. 
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2.  PREVIOUS  MEASUREMENTS 


A  survey  of  the  literature  reveals  a  vast  number  of  experimental  mea¬ 
surements  of  the  emission  spectra  of  water.  Until  recently  the  primary  Inter¬ 
est  in  the  spectra  was  the  infrared  signatures  emanating  from  rocket  exhaust 
plumes.  In  a  number  of  papers  on  the  subject,  Ludwig  and  co-workers  have  com¬ 
piled,  reworked  and  theoretically  analyzed  the  low  resolution  infrared  spectra 
of  hot  water  vapor.8,8' ^  An  example  of  their  spectra  is  shown  in  Fig.  2. 

The  spectral  region  of  interest  in  our  study  is  narrow,  944  cm-1  to 
948  cm-1,  encompassing  the  P(16),  P(18)  and  P(20)  10.6  u  CO2  laser  transition. 
Absorption  in  this  region  is  due  primarily  to  rotational  transitions  in  the 
water  molecule.  Hie  spectra  shown  in  Fig.  2  reveals  some  structure  in  this 
wavelength  region  of  interest  but  the  25  cm-1  resolution  is  not  sufficient  to 
ascertain  the  absorption  coefficient  at  a  specific  laser  transition. 

The  temperature  dependence  of  the  absorption  coefficient  can  be  esti¬ 
mated  from  low  resolution  measurements  as  shown  in  Fig.  3.  The  higher  temper¬ 
ature  values  (T  >  1000  K)  come  from  the  Ludwig  empirical  correlation. 8  The 
intermediate  values  come  from  Penner's  measurements.12'13 

In  the  lower  temperature  domain,  the  absorption  is  a  continuum.  The 
absolute  coefficient  han  a  quadratic  dependence  on  H2O  pressure.  The  absorp¬ 
tion  at  the  lower  temperature  has  been  attributed  to  water  dimers,111  (HoO)2. 
These  molecules  are  weakly  bound  and  dissociate  at  a  relatively  low  tempera¬ 
ture.  At  room  temperature  the  absorption  coefficients  have  been  measured  at 
the  lower  frequencies  of  interest  by  Peterson  and  co-workers16  and  the  absorp¬ 
tion  coefficients  are  found  to  be  very  low,  requiring  a  kilometer  of  path 
length  to  obtain  significant  absorption.  Since  water  vapor  will  condense 
into  liquid  droplets  if  the  vapor  is  saturated,  any  measurement  must  address 
the  absorption  by  liquid  water.  The  absorption  depth  of  liquid  water  at 
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10.6  p  is  only  several  microns.  If  liquid  water  is  present  in  the  gas  mix¬ 
ture  then  the  measured  absorption  coefficient  will  be  erroneously  high;  thuR, 
the  dotted  line  labeled  "Value  near  saturation"  in  Pig.  3  represents  the  maxi¬ 
mum  absorption  to  be  expected  in  pure  water  vapor. 

Prior  to  this  study  there  were  two  high  tenperature  laser  absorption 
coefficient  measurements,  one  conducted  by  Fowler  and  co-workers  at  United 
Technologies  Research  Center*  and  another  at  Physical  Sciences  sponsored  by 
NASA  Marshall  Space  Flight  Center  under  Contract  NAS8  33097. *4 

In  the  UTRC  measurements  a  water  vapor-hydrogen  gas  mixture  was  heated 
in  a  water  cooled  stainless  steel  flow  reactor  by  a  7  kW  CW  CC>2  laser.  Obser¬ 
vations  were  made  perpendicular  to  the  heating  laser  by  scanning  through  the 
hot  plasma.  T3ie  absorption  coefficients  obtained  in  these  experiments  are 
shown  in  Fig.  4,  together  with  their  predicted  values.  Hieir  measurements 
yielded  surprisingly  larqe  absorption  coefficients,  more  than  10  times  larger 
than  those  predicted  by  a  rigid  rotation  model  of  the  water  molecule.  If  ac¬ 
curate,  the  higher  value  of  the  absorption  coefficient  would  make  water  vapor 
the  prime  candidate  as  a  seed  absorber  in  a  CW  laser  heated  thruster  propul¬ 
sion  system. 

Physical  Sciences  conducted  limited  shock  tube  experiments  to  check 
the  validity  of  using  the  Ludwig  empirical  absorption  coefficients  for  water 
vapor  in  conputer  codes  developed  for  NASA  to  predict  the  performance  of  laser 
heated  thrusters.^4  The  results  of  the  experiments,  conducted  in  a  shock  tube 
using  an  untuned  CW  C02  laser,  are  shown  in  Fig.  5.  The  absorption  coeffi¬ 
cient  obtained  was  50%  lower  than  the  Ludwig  empirical  values.  Bius,  there 
existed  an  order  of  magnitude  difference  in  the  measured  values  of  the  absorp¬ 
tion  coefficient  of  water  vapor  at  10.6  microns.  Wie  present  study  seeks  to 
reconcile  these  differences. 
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3.  EXPERIMENTAL  TECHNIQUE 


3.1  Advantages  of  Shock  Tube  Technique 

T..S.  shock  tube  is  a  well  understood  experimental  tool  and  is  ideal 
for  creating  the  pressures  and  temperatures  required  for  the  water  vapor  mea¬ 
surements.  Using  a  furnace  to  produce  the  required  temperatures  becomes  very 
difficult  as  the  temperature  is  increased,  and  3000  (C  is  generally  regarded  as 
a  practical  upper  limit.  Electrical  discharges  can  reach  the  required  temper¬ 
atures  (up  to  5000  K),  but  the  resulting  gas  mixture  is  spatially  nonuniform. 
In  addition,  the  assumption  of  equilibrium  chemistry  in  the  resulting  dis¬ 
charge  is  a  poor  assumption  which  complicates  the  data  analysis.  A  flame  can 
be  used  to  create  the  required  conditions,  but  it  also  exhibits  nonuniform  gas 
mixtures.  In  fact,  flames  usually  contain  large  regions  of  cold  gas  which 
makes  the  Interpretation  of  absorption  data  difficult.  The  shock  tube  takes 
the  gas  mixture  almost  instantly  to  a  known  high  temperature  and  pressure. 

Test,  times  up  to  a  millisecond  in  duration  can  be  attained,  this  Is  more  than 
that  required  to  simulate  the  heating  process  in  a  laser-heated  thruster. 

3.2  Shock  Tube  Operating  Principles 

As  shown  in  Fig.  6,  a  shock  tube  consists  of  high  pressure  section 
and  low  pressure  sections  separated  by  a  diaphragm.  When  the  diaphragm 
bursts,  the  high  pressure  gas  flowing  into  the  low  pressure  section  acts  as  a 
piston,  compressing  and  accelerating  the  low  pressure  gas  ahead  of  it.  Com¬ 
pressing  the  low  pressure  gas  is  accomplished  by  a  shock  wave  which  travels 
ahe.'.d  of  the  piston.  When  the  shock  wave  reaches  the  end  of  the  shock  tube, 
where  the  test  section  ts.  It  reflects  off  the  end  wall  and  reshocks  the  pre¬ 
viously  shocked  gas,  compressing  more  and  stopping  it.  It  is  usual  to  specify 
the  conditions  in  the  high  pressure  region  (the  driver  gas)  by  the  subscript 
4,  conditions  in  the  undisturbed  low  pressure  section  (the  driven  gas)  by  sub¬ 
script  1,  the  gas  conditions  behind  the  incident  shock  wave  by  subscript  2 
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SHOCK  TUBE  OPERATION 


Wave  Diagram 


Pressure  Distribution  at  time  t^ 


Temperature  Profile  at  time  tj 


Fig.  6  Shock  tube  operation. 
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and  behind  the  reflected  shockwave  by  subscript  S.  Ibis  is  shown  in  the  posi¬ 
tion  versus  time  plot  (Fig.  6).  Measurements  are  made  in  the  test  section  in 
regions  labeled  2  and  5. 

The  shock  wave  heats  the  processed  (or  driven)  gas,  resulting  in  an 
enthalpy  per  unit  mass  Ah,  which  is  easily  and  accurately  calculable  when  the 
shock  velocity  Us  is  known,  Hie  approximation  Ah  ■  VB^/2  is  accurate  to  with¬ 
in  10%  for  shock  wave  Mach  numbers  greater  than  four.  Hie  enthalpy  of  the 
processed  gas  is  approximately  doubled  when  the  shock  wave  is  reflected  off 
the  end  wall.  High  tenqperatures  require  a  high  enthalpy/mass,  thus  implying 
high  shock  velocities.  Hie  velocity  of  the  gas  piston  increases  as  the  pres¬ 
sure  ratio  across  the  diaphragm  increases  but  cannot  exceed  a  value  of  a  few 
times  the  initial  sound  speed  in  the  driver  gas  (Uptaton  c  2/< Y— 1 )  a,  where  y 
is  the  specific  heat  ratio  and  a  is  the  driver  gas  sound  speed),  even  with  an 
infinite  pressure  ratio.  Hierefore,  driver  gases  with  high  sound  speeds  are 
desired.  Most  shock  tube  driver  gases  are,  therefore,  hydrogen  or  helium  and 
in  the  PSI  shock  tube,  room  temperature  hydrogen  at  up  to  200  atmospheres 
pressure  is  used. 

3. 3  Operating  Constraints 

The  initial  pressure  in  the  low  pressure  section  is  determined  by  the 
allowable  range  for  the  partial  pressure  of  water  vapor.  Hie  low  pressure 
limit  is  determined  from  the  requirements  that  any  water  vapor  on  the  shock 
tube  walls  must  not  affect  the  measured  water  vapor  mole  percent  of  the  gas  in 
the  tube.  At  low  partial  pressures  of  water  vapor  there  may  be  a  monolayer  of 
HjO  vapor  due  to  chemisorption,  and  we  have  estimated  that  van  der  Waals 
forces  might  also  yield  a  monolayer  of  water  vapor  molecules  on  the  shock  tube 
walls.  To  eliminate  any  effect  of  this  monolayer  on  the  absorption  measure¬ 
ment,  we  have  calculated  that  the  partial  pressure  of  water  vapor  must  be 
greater  than  10“3  atm. 
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Similarly,  if  the  pressure  in  the  low  pressure  section  is  raised  to 
the  point  where  the  vapor  pressure  of  water  exceeds  the  saturation  vapor  pres¬ 
sure  at  the  coldest  point  in  the  system,  water  will  condense  out  at  that 
point,  again  making  it  difficult  to  control  the  amount  of  water  in  the  tube. 

Thus,  absorption  by  water  molecules  condensed  on  the  window,  through 
which  the  absorption  measurement  will  be  made,  places  a  limit  on  the  permis¬ 
sible  operating  partial  pleasures  of  water.  Liquid  water  has  a  large  absorp¬ 
tion  coefficient  for  C&2  radiation  with  a  value  of  kv  greater  than  103  cm-1; 
therefore,  water  »wt  not  be  allowed  to  condense  onto  the  windows.  Although 
the  initial  gas  sa«i*ole  will  have  a  sufficiently  low  partial  pressure  of  water, 
the  shock  wav*»  re  is*-*  this  pressure  by  large  factors.  This  is,  of  course,  no 
problem  in  the  middle  of  the  shock  tube  where  the  gas  is  hot;  however,  the  gas 
near  the  wall,  which  is  at  this  same  high  pressure,  is  at  or  near  the  tempera¬ 
ture  of  the  wall.  It  is  important  that  the  initial  partial  pressure  of  the 
water  not  be  greater  than  the  saturation  vapor  pressure  (corresponding  to  the 
window  temperature)  at  the  wall  window  tenq>erature.  This  limit  could  have 
put  serious  constraints  on  the  operating  region.  However,  since  the  satura¬ 
tion  vapor  pressure  of  water  is  a  steep  function  of  temperature,  heating  the 
shock  tube  can  overcome  It.  Tfte  vapor  pressure  of  HjO  vapor  as  a  function 
of  temperature  is  as  follows; 

Vapor  presure  pB  =  .02  atm  @  293  K 

.12  atm  @  323  1C 
1.0  atm  @  373  K 
2.5  atm  @  400  K 
4.7  atm  §  423  K 

A  cold  boundary  layer,  if  thick  enough,  will  absorb  laser  radiation  in 
a  manner  different  from  the  uniform  hot  gas  in  the  center  of  the  tube.  It  was 
necessary  to  estimate  the  contribution  of  this  boundary  layer  to  the  overall 
absorption  measurement. 
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The  worst  case  for  boundary  layer  Interference  occurs  when  the  free 
stream  temperature  Is  600  K,  where  the  absorption  coefficient  is  a  minimum. 
Calculations  for  20%  water  in  hydrogen  at  a  pressure  of  SO  atms  indicated  a 
boundary  layer  contribution  of  12  percent.  This  contribution  gets  smaller 
rapidly  as  the  temperature  is  raised  and  is  negligible  for  temperatures 
greater  than  1000  K. 
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4.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


The  experimental  absorption  coefficient  measurements  were  performed 
behind  both  incident  and  reflected  shocks.  Hie  shock  tube  has  a  five  foot 
long  by  1.5  inch  i.d.  driver  and  a  15  foot  long  by  1.5  inch  i.d.  driven  sec¬ 
tion.  A  schematic  diagram  of  the  shock  tube  is  shown  in  Fig.  7.  Hie  optical 
measurements  are  made  through  anti-reflection  coated  zinc  sulfide  windows 
mounted  one  inch  from  the  end  wall  in  a  five  foot  long  by  1.31  inch  square 
test  section  coupled  to  the  driven  section  with  a  three  inch  long  constant 
area  transition  piece.  Shock  pressure  is  measured  by  four  piezoelectric 
transducers  located  at  one  foot  intervals  along  the  test  section.  Hie  shock 
velocity  is  determined  by  the  time  of  arrival  of  the  shock  wave  at  successive 
stations.  Hie  last  transducer  is  located  at  the  optical  port  to  allow  for  a 
direct  correlation  of  the  optical  signal  with  total  pressure. 

The  shock  heated  gas  mixtures  consisted  of  2  to  10  mole  percent  H2O  in 
either  Ar,  H2  or  N2.  Extreme  care  went  into  the  preparation  and  loading  of 
the  water  vapor/gas  mixtures  to  insure  that  the  H2O  would  not  condense  out  and 
change  the  gas  conposition.  The  tube  and  all  associated  hardware  were  main¬ 
tained  at  320  K  to  permit  the  loading  of  up  to  0.1  atm  of  H2O  vapor  without, 
condensation.  (In  actual  practice,  no  more  than  0.05  atm  of  H2<)  vapor  was 
loaded  into  the  system. )  Fresh  gas  mixtures  were  made  for  each  run.  The  sys¬ 
tem  was  evacuated  only  to  a  vacuum  of  1-10  0  to  avoid  removing  the  monolayer 
of  water  on  the  walls.  Typically,  the  system  was  exposed  to  water  vapor  by 
filling  and  evacuating  the  tube  several  times  just  prior  to  the  final  fill  be¬ 
fore  the  shock  to  further  stabilize  the  walls.  No  changes  in  pressure  with 
time  were  observed  either  after  making  a  mix  and  storing  it  in  the  mixing  tank 
or  after  loading  the  gas  in  the  shock  tube  prior  to  a  test.  Degassed,  de¬ 
ionized  water  was  stored  in  a  flask  attached  to  an  end  of  the  mixing  tank  as 
the  source  of  water  vapor.  Typical  H2O  partial  pressures  ranqed  from  16  to 
30  torr  depending  on  the  mixture  desired.  Hie  fill  pressure  was  monitored  by 
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Drlvt-r  Pressure 


Driver  Gas,  H2 


a  Validyne  pressure  transducer  with  a  0.2  torr  sensitivity.  Other  gases  were 
taken  directly  from  the  cylinder  without  further  purification. 

The  stated  minimum  purities  were:  Ar,  99.996%,  Hj,  99.995%,  N2, 
99.99%.  To  insure  good  mixing  in  the  mixing  tank,  the  diluent  gas  was  in¬ 
jected  into  the  mixing  tank  through  a  spray  bar  running  down  the  center  of  the 
tank.  Twenty  0.040"  holes  on  the  bar  at  three  inch  intervals  created  swirling 
to  promote  rapid  mixing. 

A  single-mode  line-tuned  CW  CO2  laser  was  used  as  the  source  of  the 
laser  radiation.  A  schematic  diagram  of  the  absorption  system  is  shown  in 
Fig.  8.  The  output  beam  was  split  into  two  conqponents  by  a  salt  beam  splin¬ 
ter.  About  90%  of  the  beam  was  directed  to  a  salt  diffuser  located  in  front 
of  the  entrance  window.  This  diffuser  was  employed  to  increase  the  sensiti¬ 
vity  of  the  detection  system  by  reducing  refraction  effects.  Radiation  from 
the  exit  window  was  collected  by  a  germanium  lens,  and  focused  on  an  HgCdTe 
detector  through  a  narrow  band  10.6  p  interference  filter  (FWHM  0.15  p).  The 
filter  was  used  to  block  out  water  emission  at  all  but  the  laser  wavelengths. 
No  emission  was  observed  within  the  laser  bandpass  when  the  laser  was  turned 
off  as  long  as  the  filter  was  in  front  of  the  detector.  The  basic  sensitivity 
of  the  absorption  system  is  estimated  to  be  around  0.2%  for  a  1:1  signal  to 
noise  level. 

The  remaining  10%  of  the  laser  radiation  is  sent  through  a  10  cm  path- 
length  absorption  cell  and  into  a  Scientech  laser  power  meter  to  determine  the 
laser  wavelength.  Hie  room  temperature  laser  absorption  coefficients  of  eth¬ 
ylene  have  been  accurately  measured,7  and  by  measuring  the  ethylene  absorption 
coefficients  it  is  possible  to  identify  the  laser  transition  without  using  a 
spectrometer. 
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The  absorption  coefficient  a  was  determined  using  the  measured  optical 
transmission,  I/IQ,  total  pressure  P,  path  length  1,  the  initial  mole  fraction 
of  H2O,  Xjj2o»  and  the  calculated  fractional  dissociation  at  the  equilibrium 
temperature,  fT  by  the  following  equation: 


a 


cm 


-1 


-atm 


-1 


in  (I  /I) 
o 


f  X  P 
T  H20  atm 


1 

cm 


Teagjeratures  behind  the  incident  and  reflected  shocks  were  determined 
using  the  measured  incident  shock  velocity,  the  initial  pressure  and  gas  com¬ 
positions  as  input  parameters  to  the  PSI  full  equilibrium  real  gas  shock  pro¬ 
gram.  Agreement  between  the  calculated  and  measured  shock  pressure  ratios  was 
excellent,  as  is  illustrated  in  Figs.  9  and  10. 

Experimental  data  are  collected  and  stored  in  a  computer-controlled 
CAMAC  based  data  acquisition  system.  Hie  system  Is  built  arou  ,d  two  re'*.  * 

Waveform  Digitizers.  Each  unit  can  simultaneously  sample  fou  .nalog  signals 
at  a  1  MHz  rate,  and  has  a  10— bit  32  K  word  memory  that  will  store  eight  mil¬ 
liseconds  of  data.  After  a  test,  the  data  is  transferred  to  a  PRIME  400  com¬ 
puter  for  immediate  preliminary  analysis  and  storage. 
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5.  EXPERIMENTAL  RESULTS 

The  majority  of  the  teats  were  conducted  using  argon  as  the  diluent  gas 
instead  of  hydrogen,  since  it  is  much  easier  to  drive  a  shock  into  a  heavier 
gas.  Measurements  were  made  behind  both  incident  and  reflected  shocks  encom¬ 
passing  a  temperature  range  from  600  K  to  3700  K  a  at  pressure  from  1  to  40 
atm  in  2,  5,  and  10  mole  percent  water  vapor  in  argon  gas  mixtures.  At  the 
lower  temperatures,  little  or  no  absorption  was  observed,  particularly  behind 
the  incident  shock  where  the  pressures  are  low.  Typical  low  temperature  pres¬ 
sure  and  absorption  data  is  shown  in  Figs.  11  and  12.  At  the  highest  temper¬ 
atures,  one  third  of  the  laser  radiation  was  absorbed  by  the  hot  gas  in  the 
3.31  cm  pathlength.  Typical  high  temperature  pressure  and  absorption  data  are 
shown  in  Figs.  13  and  14. 

At  several  temperatures  the  total  pressure  behind  the  shock,  and  the 
partial  pressure  of  water  in  the  gas  was  varied  sufficiently  to  assess  the  ef¬ 
fect  of  pressure  and  water  vapor  concentration  on  the  absorption  coefficient 
on  each  of  the  three  laser  transitions  of  interest.  For  example,  there  was  no 
change  in  the  value  of  the  absorption  coefficient  for  total  pressures  ranging 
from  5  to  40  atm  at  2250  K  on  any  of  the  laser  transitions  (see  Fig.  15).  At 
the  same  temperature,  plots  of  in  ID/I  vs  PH2ol  have  a  constant  slope  indica¬ 
ting  that  the  absorption  depends  solely  upon  the  line  density  of  molecules  In 
the  laser  path.  No  significant  variations  in  the  HyO  absorption  coefficient 
were  observed  among  the  three  laser  transitions  for  a  given  temperature  within 
experimental  uncertainty.  (See  Figs  16,  17,  18.)  Raw  data  collected  on  each 
of  the  laser  transitions  is  shown  in  Figs.  19,  20  and  21. 

Nine  tests  were  conducted  using  10%  H20/H2  gas  mixtures  between  600  K 
and  1400  K.  These  H20  absorption  coefficients  fall  within  H20/Ar  data  col¬ 
lected  in  the  same  temperature  range.  Two  tests  with  10%  H20/N2  gas  mixtures 
were  also  conducted  on  the  P(18)  C02  transitions  and  these  also  overlap  the 
H20/Ar  data.  All  the  data  collected  along  with  the  Ludwig  empirical  curve  fit 
i s  plotted  in  Fig.  22. 
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Typical  low  temperature  pressure  trac 


Fig.  13  Typical  high  temperature  pressure  trace 
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Fig.  21  Absorption  coefficient  vs.  temperature  for  (20)  laser  line 


6.  DISCUSSION 


The  experimental  observation  that  changes  In  total  gae  pressure,  the 
partial  pressure  of  the  water,  or  the  collision  partner,  along  with  small 
changes  In  the  laser  wavelength,  do  not  significantly  change  the  value  of  the 
absorption  coefficient  of  H2O  vapor/gas  mixtures  over  the  ranges  of  pressures 
and  temperatures  encotq>assed  by  these  measurements  is  not  surprising  for  the 
following  reason. 

A  continuum  absorption  will  be  obtained  at  any  temperature  if  the 
lines  involved  in  the  absorption  are  broadened  into  each  other.  As  the  tem¬ 
perature  rises  the  density  of  states  increases  so  that  the  degree  of  broaden¬ 
ing  required  to  complete  the  overlap  decreases. 

The  positions  of  the  high  temperature  water  lines  are  not  known  with 
sufficient  accuracy  to  allow  a  line  by  line  calculation  of  the  absorption  co¬ 
efficient.  Ludwig  and  his  coworkers  have  conducted  many  measurements  of  water 
vapor  emission  from  plumes  and  flames,8  empirically  correlated  all  available 
low  resolution  spectra, c>' and  obtain  an  expression  for  the  temperature  and 
pressure  dependence  of  the  absorption  coefficient  averaged  over  a  25  cm-1 
spectral  range. 

The  band  model  used  by  Ludwig  determines  an  average  rotational  line 
spacing  vs.  temperature  as  shown  in  Fig.  23.  His  equation  for  the  average 
spacing  as  a  function  of  temperature  is  d(T)  cm-1  »  exp(-0. 00106  TK  +  1.21). 
The  primary  broadening  mechanism  under  the  experimental  conditions  is  colli- 
slonal  broadening.  The  collisional  broadening  parameters  for  all  the  species 
of  interest  are  known  at  room  temperature, 1 1  and  using  these  values  it  is  a 
simple  matter  to  approximate  the  pressure  required  to  broaden  a  line  to  an  ex¬ 
tent  that  the  linewidth  vt  equals  the  average  line  spacing  d.  Above  this  pres¬ 
sure  the  absorption  will  be  a  continuum.  These  plots  for  the  gases  tested  in 
the  study  are  shown  sb  Fig.  24. 
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For  the  pressure  conditions  required  in  the  operation  of  a  Cw  laser- 
heated  thruster,  this  simple  approach  predicts  sufficient  broadeninq  to  assume 
a  continuum  absorption  mechanism.  Our  data  closely  approximates  a  curve  60% 
lower  than  Ludwig's  empirical  fit.  Ibis  is  not  unexpected  either  since  there 
is  some  structure  in  the  high  temperature  water  spectra  in  the  10.6  micron 
region.  Low  resolution  spectra  recorded  in  the  944-948  cm~^  region  of  the 
spectra  appear  to  be  slightly  below  Ludwig's  averaged  value;  however  the  re¬ 
spective  high  resolution  spectra  could  have  varied  in  either  direction. 

No  absorption  was  observed  below  1000  K.  This  results  from  a  combina¬ 
tion  of  short  path  length,  low  pressure,  and  low  absorption  coefficients.  In 
several  other  measurements  it  has  been  observed  that  the  absorption  coeffi¬ 
cient  increases  quadratically  with  increasing  pressure.  Ibis  continuum-like 
phenomenon  has  been  attributed  to  absorption  by  water  dimers.  The  weakly 
bound  complexes  dissociate  as  the  temperature  increases,  and  their  dissocia¬ 
tion  is  suppressed  by  increasing  pressure.  The  functional  form  of  the  low 
temperature  absorption  has  been  described  by  Roberta  and  coworkers. ^ 
Combining  our  hi^i  temperature  experimental  measurements  fitted  to  a  value  60% 
of  Ludwig's  predictions,  and  Roberts  1<*»  teaq>erature  observations,  a  plot  of 
the  temperature  dependence  of  the  absorption  coefficient  of  H20  can  be  gener¬ 
ated  for  all  possible  operating  regimes  of  a  CW  laser-heated  thruster.  Ibis 
plot  is  shown  in  Fig.  25. 
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Calculated  collision  broadening  boundaries  for  test  aas  mjvfuro 
versus  temperature. 
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Fig.  25  psi  reconmended  values  for  HO  absorption  coefficient. 
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7.  WATER  AS  AM  ABSORBER  IN  THE  CW  ROCKET 


The  low  value  of  the  absorption  coefficient  of  water  below  1500  K  in¬ 
dicates  that  a  laser  heated  thruster  using  a  simple  h2o/h2  propellant  system 
may  have  a  startup  problem,  ms  will  not  be  known  until  it  is  possible  to 
design  and  test  a  laboratory  scale  thruster.  The  reason  is  that  the  presence 
of  liquid  water  droplets  in  the  gas  stream  will  cause  a  significant  increase 
in  the  absorption  coefficient,  and  it  might  be  possible  to  heat  the  propellant 
in  a  rather  short  absorption  chamber  by  superheating  the  droplets  and  genera¬ 
ting  the  steam  at  a  temperature  where  the  gas  phase  water  absorbs  strongly. 

If  we  neglect  the  potential  start  up  problem,  we  can  investigate  the 
performance  of  a  prototype  laser  heated  thruster  using  the  PSI  laser  heated 
thruster  codes. ^  Figure  26  indicates  the  specific  Impulse  and  absorption 
length  for  a  variety  of  H20/H2  propellant  ratios.  The  diagonal  family  of 
curves  represent  the  absorption  lengths  required  for  various  propellant  mix¬ 
tures  vs  chamber  pressure.  The  gas  is  assumed  to  be  either  at  3500  or  4000  K. 
The  more  horizontal  family  of  curves  represents  the  specific  impulse  obtained 
under  these  conditions.  This  figure  indicates  that  a  reasonably  compact 
thruster  operating  at  3500  to  4000  K  and  above  10  atmospheres  pressure  can 
yield  specific  impulses  around  1000  s. 

The  laser  power  required  and  absorption  chamber  heating  loads  for  a 
small  2  cm  radius  laser  heated  thruster  are  shown  in  Fig.  27.  This  figure 
shows  that  with  a  30-50  kW  CW  C02  laser  system,  this  size  thruster  will  pro¬ 
duce  1  pound  of  thrust  with  an  average  heating  load  of  only  40  to  80  W/cm2. 
This  is  not  an  unreasonable  power  level  for  state  of  the  art  C02  lasers,  thus 
demonstration  experiments  for  a  1  pound  thruster  can  in  principal  be  con¬ 
ducted. 
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8.  AMMONIA  AS  A  PROPELLANT  AND  ABSORBER 

As  stated  in  the  previous  section,  there  might  be  a  start  up  problem 
in  a  laser  heated  thruster  using  an  H2O/H2  propellant  system.  In  anticipation 
of  this  possibility,  we  have  reviewed  the  available  experimental  data  for  am¬ 
monia,  a  molecule  which  strongly  absorbs  CO2  laser  radiation  at  room  tempera¬ 
ture. 


Room  temperature  laser  absorption  coefficient  measurements  of  ammonia 
have  been  conducted  for  the  laser  transitions  of  interest7  .  The  absorption 
coefficients  ares  aP(16)  *  0.12,  aP(18)  ■  0.14,  and  aP(20)  »  0.49  cm-1  atm-1. 
The  only  high  temperature  laser  absorption  measurements  found  are  those  con¬ 
ducted  at  UTRC*  and  their  values  are  used  in  the  following  discussion. 

It  is  desirable  to  use  NH3  as  the  propellant  in  a  CD  laser  heated 
thruster  because  NH3  is  easily  stored  as  a  liquid  and  leads  to  a  product  (N2  + 
SHj)  with  a  low  average  molecular  weight.  If  the  NH3  remains  intact  at  higher 
temperatures,  it  could  serve  as  an  effective  lew  temperature  absorber.  How¬ 
ever,  an  equilibrium  calculation  of  the  mole  fraction  of  NH3  resulting  from  an 
intially  pure  NH3  gas,  shown  in  Pig.  28,  indicates  that  from  an  equilibrium 
standpoint,  NH3  does  not  exist  in  appreciable  quantities  at  temperatures  above 
1000  K.  Therefore,  only  if  the  heating  process  remains  nonequilibrium,  will 
the  NH3  remain  intact.  The  kinetics  for  the  decomposition  of  NH3  are  reason¬ 
ably  well  known  and  are  fairly  slow  at  temperatures  below  3000  K. 

We  have  performed  several  calculations  to  determine  the  rate  of  decom¬ 
position  of  NH3  at  various  temperatures  and  pressures  in  order  to  make  a  pre¬ 
liminary  assessment  of  its  use  as  a  low  tenqperature  absorber.  If  NH3  could 
exist  in  significant  quantities  to  tenqperatures  where  absorption  by  H2O  could 
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take  over,  NH3  would  indeed  be  a  useful  propellant  additive.  In  Fig.  29  we 
show  the  results  of  these  calculations  where  we  have  calculated  the  time  his¬ 
tory  of  the  NH3  concentration  after  sudden  heating  (such  as  in  a  shock  tube) 
to  various  temperatures  at  various  pressures. 

The  solid  curve  shows  the  NH3  history  at  10  atmospheres  and  3000  K. 

In  these  calculations  we  have  made  the  area  of  the  flow  constant;  the  gas 
cools  rapidly  to  2500  K  due  to  the  Initial  breakup  of  the  NH3  to  form  radi¬ 
cals.  Hie  rate  of  removal  of  NH3  after  about  1  us  is  controlled  by  the  reac¬ 
tion 


N  +  NH  *  N2  +  H 

which  is  limited  by  the  low  concentration  of  nitrogen  atoms.  At  higher  pres¬ 
sure  (40  atm)  the  rate  of  removal  is  even  slower  as  seen  from  the  dashed  curve 
in  Fig.  29.  This  is  because  the  nitrogen  atom  concentration  is  in  quasi- 
steady  state  with  N2  and  is  thus  reduced  by  a  factor  of  the  square  root  of  the 
pressure.  With  the  addition  of  10%  water,  the  rate  of  NH3  removal  is  in¬ 
creased  by  a  factor  of  two  as  is  Shown  by  the  dash-double-dot  curve.  This  is 
due  to  the  introduction  of  a  new  mechanism 

N  +  OH  ♦  NO  +  H 
O  +  NH  ♦  NO  +  H 
N  +  NO  ♦  N2  +  0 

This  mechanism  still  depends  upon  the  presence  of  nitrogen  atoms  so  that  it 
will  scale  in  the  same  way  as  the  decomposition  without  H20. 

At  2000  K  the  decon$>osition  rate  is  very  slow.  Ibis  is  Indicated  in 
Fig.  29  by  the  circle  at  10  ms.  At  4000  K  the  NH3  decomposition  rate  is  fast 
as  la  Indicated  by  the  dash-dot  curve.  From  these  curves  we  can  expect  that 
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for  a  device  with  a  characteristic  heating  time  of  about  1  ms,  the  absorption 
due  to  NH3  would  drop  off  slowly  below  2500  K  and  rapidly  above  3000  K.  This 
is  qualitatively  what  was  seen  in  the  ejqperiments  of  Fowler,  et  al,4  at  UTRC. 
The  results  of  their  e>q>erimenta  are  reproduced  in  Fig.  31.  Hie  rise  in  the 
absorption  at  temperatures  above  1000  K  is  not  understood  but  could  be  due  to 
the  increase  in  the  population  of  absorbing  hot  bands  of  NH3  or  radicals.  Hie 
fall  is  undoubtedly  due  to  equilibration  of  the  conposition. 

With  the  kinetics  of  NH3  dissociation  established  in  Fig.  29,  we  would 
like  to  compare  the  dissociation  time  to  the  flow  time  in  the  heating  region 
of  a  laser-heated  thruster  to  determine  if  NH3  would  make  an  effective  ab¬ 
sorber  when  used  in  a  thruster.  An  estimate  of  the  flow  time  at  a  temperature 
level  in  a  laser-heated  thruster  can  be  made  by  combining  the  flow  velocity 
expected  with  the  absorption  length. 

Suppose  we  have  a  thruster  operating  to  produce  a  thrust  Tb  at  a  given 
specific  impulse  Iap.  Then  the  mass  flow  rate  is  m  -  Th/g  Isp  where  g  is  the 
gravitational  constant.  If  A  is  the  thruster  croaa  sectional  area  and  p  the 
density  at  the  cross  section  of  interest,  then  the  velocity  there  is 

u  -  Th/(g  I8p  pA) 

The  absorption  coefficient  may  be  expressed  as  an  absorption  cross-section  o 
times  the  number  density  of  NH3,  call  it  n.  Hius 

k  ■  on  ■  crp/m 

where  m  is  the  mass  of  an  ammonia  molecule.  By  combining  u  and  k  we  obtain 
the  flow  time 
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t  .  1_  ■  15_ 

uk  a 


where  R  is  th«  cross-section  radius,  this  tins  is  assn  to  be  independent  of 
pressure,  but  depends  on  temperature  through  the  absorption  cross-section  a. 

It  is  directly  proportional  to  area  (r2),  and  Isp,  and  inversely  proportional 
to  thrust,  Th.  High  thrust  decreases  the  time  and  large  area  or  Isp  Increases 
it. 


To  make  an  estimate  for  NH3  we  need  to  find  o.  For  this,  we  use  Fig. 
30  of  Ref.  4,  reproduced  as  Fig.  31  here.  This  is  a  plot  of  the  calculated 
absorption  coefficient  of  1  atm  of  NH3  in  10  atm  of  The  upper  curve  is 

for  a  frosen  mixture,  and  so  Is  only  sensitive  to  the  pressure  of  NH3.  The 
following  table  gives  the  values  of  k  from  the  upper  curve  of  Fig.  31,  the 
number  density  of  1  atm  of  NH3,  and  their  ratio  a  ■  k/n,  against  tenqperature. 


T(103K) 

k ( cm" 1 ) 

n(cm" 

h 

0.3 

0.61 

2.45 

X 

1019 

2.49 

X 

10-2° 

1 

0.35 

7.34 

X 

to1® 

4.77 

X 

10-20 

2 

0.Q70 

3.67 

X 

103® 

1.91 

X 

10-20 

2.5 

0.032 

2.94 

X 

to18 

1.09 

X 

10-20 

3 

0.014 

2.45 

X 

10™ 

5.72 

X 

10-21 

4 

0.0037 

1.84 

X 

lo’8 

2.02 

X 

10-21 

5 

0.0011 

1.47 

X 

10’8 

7.49 

X 

10-22 

6 

0.00035 

1.22 

X 

1018 

2.78 

X 

10-22 

The  mass  of  an  ammonia  molecule  is  2.85  x  10“23  g  so  the  expression 
for  t  becomes 
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ABSORPTION  PER  CM  FOR  1  ATM  NH3 


1000  2000  3000  4000  5000  6000 

T-°K 


Fig.  31  UTRC  calculated  NH,  absorption  coefficient 


t  - 


8.8  x  TQ-20r2  *sp 
o  Th 

where  o  is  in  cm2,  r  is  in  cm  and  ’Hi  is  in  dynes,  while  t  and  I8p  are  in 
seconds,  'ftiis  can  be  used  to  estimate  the  flow  time  in  a  thruster  of  radius 
R,  thrust  Th,  and  specific  impulse  ISp*  at  the  temperature  where  the  absorp¬ 
tion  cross-section  for  ammonia  is  o. 

As  an  example,  consider  a  thruster  developing  100  pounds  thrust 
(4.45  x  1  (P  dynes)  at  a  specific  impulse  of  1000  s.  Then 

t  »  2  x  10-24  R2/o 

Using  the  values  for  o  from  the  above  table,  we  find 


R-1 

era 

R=3  cm 

T( 1 03K) 

t/R2(s/cm2) 

t(s) 

t(s) 

1 

4.2 

X 

in 

1 

o 

e— 

4.2  x 

10-5 

3.8  x  10"4 

2 

1.0 

X 

10-4 

1.0  x 

O 

1 

9.0  x  10“4 

2.5 

1.8 

X 

IQ’4 

1.8  x 

o 

1 

As 

1.8  x  10-3 

3 

3.5 

X 

1 

© 

3.5  x 

* 

I 

o 

r- 

3.2  x  10-3 

4 

9.9 

X 

IQ-4 

9.9  x 

10-4 

8.9  x  IQ-3 

At  3000  K  Fig.  30  and  the  above  table  indicate  that  the  NHj  would  be 
more  than  60%  dissociated  for  this  thruster  if  it  had  a  radius  of  1  cm,  while 
it  would  be  90%  dissociated  if  it  had  a  radius  of  3  cm.  However,  at  2500  K, 
it  is  only  20%  dissociated  for  a  radius  of  1  cm  and  30%  dissociated  for  a 
radius  of  3  cm.  So  it  appears  that  the  kinetics  of  NH3  decomposition  is  slow 


enough  that  NHj  will  remain  intact  and  be  an  effective  low  temperature  absor¬ 
ber  up  to  2500  K. 

Based  upon  the  above  discussion  we  conclude  that  ammonia  in  a  laser 
powered  thruster  might  serve  as  an  effective  absorber  during  the  early  heating 
phases.  Since  water  vapor  is  not  a  good  absorber  at  low  temperatures 
(<1500  K)  and  ammonia  is.  a  NH3/H2O/H2  mixture  is  attractive  as  a  10.6  pm  ra¬ 
diative  absorber  that  operates  from  room  temperature  to  4000  K,  where  water 
vapor  is  the  absorber  above  2500  K. 
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9.  SUMMARY  AND  CONCLUSIONS 

No  significant  variation  in  the  H20  absorption  coefficient  was  ob¬ 
served  as  a  function  of  laser  wavelength,  water  concentration,  total  pressure, 
or  collision  partner  within  the  narrow  spectral  range  from  944  to  94H  cm-1 
covered  in  the  measurements.  Wiese  observations  suggest  that  the  water  lines 
are  sufficiently  broadened  to  act  as  a  continuum  absorber  under  conditions  to 
be  found  in  a  laser-heated  rocket  thruster,  Wie  laser-measured  high  temper¬ 
ature  absorption  coefficients  are  50  percent  lower  than  the  value  obtained 
from  the  Ludwig  empirical  curve  fit  to  low  resolution  data. 

For  a  practical  CW  laser-heated  thruster  it  is  estimated  that  a  pro¬ 
pellant  absorption  coefficient  of  10“^  atm“^-cm”^  will  be  required. ^»4« ^ 

The  absorption  coefficient  of  water  vapor  does  not  reach  this  value  until  the 
gas  temperature  exceeds  1500-1700  K.  High  temperature  operation  will  not  pre¬ 
sent  problems,  but  there  could  be  a  potential  start-up  problem  for  a  simple 
H20/H2  propellant  system,  and  other  low  temperature  absorbers  may  be 
required. 

Ammonia  absorbs  C02  laser  radiation  strongly  from  room  temperature  to 
around  2500  K.  Although  equilibrium  calculations  predict  substantial  dissoci¬ 
ation  at  elevated  temperatures,  the  kinetics  suggests  that  the  decomposition 
of  ammonia  below  2500  K  is  sufficiently  slow  as  to  make  NH3  a  very  attractive 
low  temperature  absorber.  Addition  of  NH3  will  eliminate  any  low  temperature 
start  up  problems,  A  tertiary  propellant  mixture,  NH3/H20/H2,  therefore 
should  absorb  C02  laser  radiation  effectively  from  room  temperature  to  above 
4000  K. 
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